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Abstract 

The site percolation problem is studied on d-dimensional generalisations of the Kagome 
lattice. These lattices are isotropic and have the same coordination number q as the hyper- 
cubic lattices in d dimensions, namely q = 2d. The site percolation thresholds are calculated 
numerically for d — 3, 4, 5, and 6. The scaling of these thresholds as a function of dimension d, 
or alternatively q, is different than for hypercubic lattices: p c ~ 2/q instead of p c ~ l/(q— 1). 
The latter is the Bethe approximation, which is usually assumed to hold for all lattices in high 
dimensions. A series expansion is calculated, in order to understand the different behaviour 
of the Kagome lattice. The return probability of a random walker on these lattices is also 
shown to scale as 2/q. For bond percolation on d-dimensional diamond lattices these results 
imply Pc ~ 1/(5- 1). 

1 Introduction 

The Kagome lattice is one of the most interesting lattices in two dimensions. It is one of the 
eleven Archimedean tiling lattices, where all the vertices are of the same type (see e.g. Weisstein 
1997). In the case of the Kagome lattice each vertex touches a triangle, hexagon, triangle, and 
a hexagon. All these polygons are regular. Moreover, the Kagome lattice is closely related to 
the other lattices in two dimensions. The sites of the Kagome lattice correspond to the bonds 
of the honeycomb lattice, which in turn is the dual of the triangular lattice. Therefore, since 
the bond percolation threshold of the honeycomb lattice is 1 — 2sin(7r/18), the site percolation 
threshold of the Kagome lattice is 1 — 2sin(7r/18) = 0.6527036. . . too (Sykes and Essam 1964). 
The bond percolation threshold is not known exactly, but has been calculated numerically with 
high precision to be 0.524 405 3 ± 0.000 000 5 (Ziff and Suding 1997). 

Although these percolation thresholds have been known for quite some time, it is not clear 
why this site percolation threshold is high, compared to other lattices. For instance, it is much 
higher than the threshold 0.592 746 ± 0.000 000 5 for the square lattice (Ziff and Sapoval 1986, 
Ziff 1992), although its coordination number q = 4 is equal to that of the Kagome lattice. What 
is more striking is that even the pentagonal lattice, which has a low average coordination number 
of q = 3|, has a lower site percolation threshold 0.6471 ± 0.0006 than the Kagome lattice (Van 
der Marck 1997). In other words, the site percolation thresholds are not ordered according to 
the coordination number q. This runs contrary to common intuition, which leads one to expect 
that a lattice with a higher connectivity has a lower percolation threshold. If one searches for 
general formulas that correlate percolation thresholds with dimension and coordination number, 
the Kagome lattice therefore poses a problem. Galam and Mauger (1996) introduced different 
classes of lattices to avoid this problem, and they used different correlations for these classes. 
Although this enabled them to derive good correlations, it prompts the question why certain 
lattices belong to one class and others to another class. 







An analogue in three dimensions was found recently: there is a lattice with q = 6 and perco- 
lation threshold 0.3898 ± 0.0008 (Van der Marck 1997b). Compared to the simple cubic lattice, 
which also has q = 6, but a threshold of 0.31 14 ± 0.0004, this threshold is much higher. It is even 
higher than the threshold for several lattices with coordination number q — 5 (Van der Marck 
1997). 

In this paper, a generalisation of the Kagome lattice to higher dimensions is given (Sec. |^), 
and numerical calculations of the site percolation thresholds for 3, 4, 5, and 6 dimensions are 
presented (Sec. ||). The scaling of these thresholds as a function of dimension appears to be 
different than for hypercubic lattices. For the latter, Gaunt, Sykes and Ruskin (1976) calculated 
a series expansion in l/(2d — 1), where d is the number of dimensions. The leading term in their 
series is p c (d) = l/(2d— 1), which is the so-called Bethe approximation. This approximation 
holds exactly for Bathe lattices (see e.g. Stauffer and Aharony 1992). For hypercubic lattices, 
the approximation underestimates the percolation threshold in low dimensions, but improves in 
accuracy in higher dimensions. In Sec. || the series expansion for the <i-dimensional Kagome lattices 
is studied. It is suggested that in this case the leading term in the series is 1/d, not l/(2d — 1). 

The return probability of a random walker on d-dimensional Kagome lattices is discussed in 
Sec. ||. Ishioka and Koiwa (1978) conjectured that this probability is a good estimator for the 
percolation threshold. Indeed it is shown here that this return probability also scales as 1/d. 
Section || contains a discussion of the results. These results give some insight into the problem 
why the 2-dimensional Kagome lattice has a high site percolation threshold. 

2 Kagome lattices in d dimensions 

The Kagome lattice can be defined in d dimensions as follows. The lattice has a (d + l)-point 
basis, and these 

points form a regular d-dimensional polytope. All the points of this basis are direct neighbours 
of each other. Let us denote the basis points by b^, for i = 0, . . . , d. The lattice can be constructed 
by translation of the basis in ^d(d + 1) directions. These translations are given by the vectors 
2(bj — hi) for j 7^ i. This is a dependent set of vectors. One can select a minimal set of d vectors 
by setting e.g. i = 0, j = 1, . . . , d. 

A site of the lattice can be identified by its number i within the base polytope (runs from 
to d), and the translation x with respect to a reference position. Consider a site {i,x}. This 
site has d neighbours in the same polytope, {j ^ *,x}, and another d neighbours in adjacent 
polytopes, {j z,x + 2(b^ — bj)}. Therefore the coordination number of this lattice is q = 2d. 
Also, because there is no preferential direction in the construction, all the directions are equivalent 
for this lattice, i.e. it is an isotropic lattice. 

So the Kagome lattice resembles the cubic lattice in the sense that both are isotropic d- 
dimcnsional lattices with coordination number q = 2d. However, we already know that in two and 
three dimensions, the site percolation thresholds of these lattices are distinctly different. 

In two dimensions the site percolation problem on the Kagome lattice is equivalent to the bond 
percolation problem on the honeycomb lattice (Sykes and Essam 1964). This was shown by means 
of the star-triangle transformation. Analogous to the star-triangle transformation, one can use a 
'star-tetrahedron' transformation in three dimensions, see Fig. |^. The site percolation problem 
on the tetrahedron ABCD, i.e. the three-dimensional Kagome lattice, is equivalent to the bond 
percolation problem on the dashed lattice. This dashed lattice is the diamond lattice, for which 
Pcb — 0.3893 ± 0.0003 has been calculated (Van der Marck 1997). Generalising to d dimensions, 
the site percolation on the d-dimensional Kagome lattice is equivalent to bond percolation on the 
d-dimensional diamond lattice, which has q = d + 1. 

The Kagome lattice defines a tiling of d-dimensional space. In two dimensions, the Kagome 
lattice defines a tiling of the plane in terms of a regular triangle and a regular hexagon. One 
can construct this hexagon from the base triangle: take three neighbouring base triangles, as in 
Fig [l]b. These three form a larger regular triangle. The hexagon appears when one truncates 
the larger triangle by taking away the three smaller triangles. In two dimensions this is a rather 
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complicated description of the tiling, but the advantage is that one can use an identical procedure 
in d dimensions. One can start with the regular base polytope of (d+ 1) points in d dimensions, 
(for i = 0, . . . , d). This polytope is shifted in d directions by the vectors 2(b^ — bo) for i = 1, . . . , d. 
This defines a larger regular polytope, which is then truncated by taking away the (d + 1) small 
polytopes. In three dimensions, e.g., the resulting polyhedron is the truncated tetrahedron (see 
e.g. Weisstein 1997). The d-dimensional space is filled with the base regular polytope and the 
larger truncated polytope. 

In view of the considerations in sections ^ and it is helpful to discuss one more important 
property of the d-dimensional Kagome lattice: two adjacent sites have d — 1 common neighbours. 
This is demonstrated by the following arguments. When the two sites are within the same base 
polytope, characterised by {i,x} and {j, x} (j ^ i), all the other d — 1 sites of that polytope 
are common neighbours. When the two sites are not within the same base polytope, they are 
characterised by {i, x} and {j, x + 2(bj — hj)}. Note that j cannot be equal to i for two adjacent 
sites. The d — 1 common neighbours of these two sites are {fc, x + 2(b^ — b^)} with the restriction 
that k ^ This is because they are one step b; — bfe away from {i,x}, and one step b^ — hj 
from {j, x + 2(hi — hj)}. 

3 Percolation thresholds 

The site percolation thresholds of the Kagome lattice in 3, 4, 5, and 6 dimensions, are show in 
Table | as a function of the lattice size. I used the method given by Stauffer and Aharony (1992). 
This method is simple to program, but in its simplest form it has the drawback that one needs to 
have all N — (d+l)L d sites of a network resident in memory. Therefore the calculation is restricted 
to relatively low values of L, the linear size of the lattice. Especially in higher dimensions this 
is rather restrictive: the highest L-value used was L — 14 for the 6-dimensional Kagome lattice. 
One can check in a scaling plot, Fig. [| that the used lattices are indeed large enough to be in the 
scaling regime, and hence allow extrapolation to infinite lattice sizes. 

For each of the lattices, percolation in only one direction, say the xi-direction, was checked. 
The lattice was defined to be percolating whenever the two boundaries {i = 0,xi = 1} and 
{i = d, xi = L} were connected. The boundary conditions in the other directions were not 
periodic. 

The percolation thresholds of a lattice of size N obey the scaling relation 

\p c (N) - Pc (oo)\ ~ ATVM). (1) 

Here the critical exponent v is 0.88 in three dimensions and 0.68, 0.57, 0.5 in 4, 5, and 6 dimensions 
respectively (Stauffer and Aharony 1992). The results quoted in the row marked oo in table ^ are 
fits of p c (N) to this scaling relation, using the last three points in the table. 

The data and the fits are shown in Fig. || for d — 3 and d = 6. From the figure it can be 
concluded that the input values for the fits are indeed within the scaling regime, and that the fits 
are therefore sensible. For d = 6 the difference p c (N) — p c (oo) is larger than for lower dimensions, 
because the number of lattice sites that are near to a boundary increases with dimension. 

Also shown in Table |l| are results for the percolation thresholds of the cubic lattice. These 
numbers agree with the literature (Stauffer and Aharony 1992) and can be seen as a check that 
the programs used in the present work are correct. 

The difference between the percolation threshold for the Kagome lattices and the (hyper-) 
cubic ones increases with dimension. This is illustrated in Fig. |]. In the next section a series 
expansion is studied, in order to understand this behaviour. 

4 Series expansion 

One could conjecture, based on the numerical estimates of the percolation thresholds, that the 
scaling of these thresholds as a function of dimension is different for the Kagome lattice than for 
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the cubic lattice. The scaling behaviour of the cubic lattice was studied by Gaunt et al. (1976). 
Their result for the percolation threshold was 

1 A 3/2 15/4 83/4 \ 

The leading term is the percolation threshold of a Bethe lattice with coordination number q = 2d. 
Equation || agrees well with the known numerical estimates for the cubic lattices for d > 3, as is 
clear from Fig. [|. 

Gaunt et al. (1976) used cluster counting to calculate a series expansion for the mean size S of 
clusters at low probabilities p. Since the mean cluster size diverges at the percolation threshold, 
they studied the radius of convergence of the series expansion of S. They argued that this radius of 
convergence will be determined predominantly by the singularity at p = p c , which enabled them to 
derive Eq. (0) . In this section these techniques are used to study the percolation thresholds of the 
Kagome lattice in higher dimensions. At the basis of the approach are the 'perimeter polynomials', 
as described by Sykes and Glen (1976). Denoting the mean number per lattice site of clusters of 
s sites by n s , we have e.g. 

«i = v i 2d > 

n 2 = dp 2 q M ~ X . (3) 

Here p is the probability that a site is occupied, and q = 1 — p the probability that a site is empty. 
One can interpret p(s) — sn s as the probability that a site is occupied by a cluster of size s. When 
one sums p(s) over all cluster sizes, one gets the probability that a site is occupied by any cluster, 
i.e. p: 

oo oo 

^2sh s = ^2p(s)=p. (4) 

s=l s=l 

The power 3d — 1 of q in the expression for n 2 indicates that all clusters of 2 sites are surrounded 
by 3d— 1 neighbours. Already we see a difference with the cubic lattice, which has n 2 — dp 2 q Ad ~ 2 . 
So even though a site on the Kagome lattice has as many neighbours as a site on the cubic lattice 
(namely 2d), a 2-site cluster has less neighbours on a Kagome lattice than on a cubic lattice. This 
is because on the Kagome lattice two adjacent sites have d — 1 common neighbours, whereas on 
the cubic lattice they have none. 

A relatively simple counting procedure by hand reveals the next few perimeter polynomials to 

be 

n 3 = ±d(4d-l)p 3 q id -\ (5) 
n 4 = ±d(5d-l){5d-2)p i q 5d - 3 , 

n 5 = [±d(6d-l)(6d-2)(6d-3)-2d(d-l)]p 5 q 6 ^ 
+2d{d~ 1) p 5 q 6d - 5 . 

In Figs. |^ and |^ the contributing graphs are depicted in two dimensions. In the expressions 
for fix, . . . ,ri4 a clear pattern seems to emerge. The coefficient of n s has factors (s + l)d — 1, 
(s + l)d — 2, . . ., and the power of q is always (s + l)d — (s — 1). 

However, this trend is broken in 715. Although the sum of the two coefficients appearing here 
still has the factors (s + l)d — 1, etc., there are now terms with different powers of q. In other 
words, in ns we encounter for the first time that there are clusters with different numbers of 
neighbours; clusters with 6d — 4 neighbours and clusters with 6d — 5 neighbours. The latter ones 
are depicted in Fig. ||. They are the type of cluster that 'bites itself in the tail'. At first sight, one 
would think that this would occur already for clusters of two sites, but in those cases the number 
of neighbours is not affected. This number is only affected for the clusters depicted in Fig. |[ On 
the cubic lattice this occurs for the first time for clusters of 3 sites. Another difference is that for 
the cubic lattice the leading power of q is 2sd — 2(s — 1), compared to (s + l)d — (s — 1) for the 
Kagome lattice. 
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The identity (Q) can be exploited to determine one more mean cluster number, namely ng, 
albeit in the limit of q — > 1. Following Sykes and Glen (1976) I substitute Eqs. (||) and (JsJ) into 
the identity (Q) , and set q — 1 — p. This yields a power series in p, of which the coefficient of 
p is 1, and the coefficients of p 2 , . . . ,p 5 vanish identically. The higher order coefficients do not 
vanish, because the number of clusters of six or more sites were not included. Because uq is the 
only missing term that can contribute terms of order p 6 , we can conclude 



n e\ q . 



±-d(7d - l)(7d - 2)(7d - 3)(7d - 4) - ^d(d - 1) 
360 3 



A (6) 



The mean size of clusters at low probabilities, S, is defined as 

3 = 



I>p( S )_l£ s2 _ s _ (7) 



E s pO) p 

Again Eqs. (||) and (||) can be used, and q = 1 — p, to derive a power series 



S(p) = 1 + 2(dp) + 2(rfp) 2 + 2(dp) 3 + 2(dp) 



2(dpf 



10 / 1 
1 + ^1-^ 



d 3 V d 



0( P b ). (8) 



The first few terms in this low density expansion of S(p) are remarkably simple. It is almost a 
geometric series, until the fifth power. The extra terms in the coefficient of (dp) 5 are due to the 
clusters that 'bite in their own tail'. 

Although the resemblance of S(p) to a geometric series is not exact, it does suggest that 1/d is 
the obvious candidate for expansion parameter. This in contrast to the cubic lattice, where it is 
l/(2d- 1), see Eq. (|). In fact, the resemblance to a geometric series suggests that a singularity 
of S(p) should occur in the vicinity of dp = 1. Therefore one can expect the percolation threshold 
to scale as 

Pc 4 (9) 

In Fig. [I| the relation p c = 1/d is shown with a dashed line. The percolation thresholds for d = 5 
and 6 are already reasonably well approximated by this relation. In an attempt to calculate the 
percolation threshold for d = 8, I computed 0.1059 ± 0.0005 for L = 5, 0.1086 ± 0.0006 for L = 6, 
and 0.1117 ± 0.0005 for L = 7. Based on the last two points one can fit the percolation threshold 
to be 0.120 ± 0.003, which is close to the value 0.125 one would expect on the basis of Eq. (0). 
However, the point for L = 5 does not lie on the same fit, indicating that these lattice sizes are not 
yet large enough. Therefore there is probably also a small systematic error in the determination 
of the value 0.120. 

When one wants to refine the scaling behaviour given by Eq. (^) with terms of the order of 1 /d 2 
and further, analogous to Eq. (|^), more terms in the series expansion would be required, plus a 
careful mathematical analysis of the radius of convergence of the series. This is beyond the scope 
of the present paper. 



5 Random walks 

It is interesting to study the return probability P r of a random walker on a d-dimensional Kagome 
lattice. Ishioka and Koiwa (1978) suggested that P r is an upper bound for the site percolation 
threshold on any lattice: P r > p c . For Bethe lattices, the equal sign holds (Hughes and Sahimi 
1982). The arguments given by Ishioka and Koiwa to support their conjecture are not exact, but 
the relation appears to work fairly well. P r and p c lie closer together, the more connected a lattice 
is. Also Sahimi et al. (1983) studied the relation between a random walker (not self-avoiding) and 
percolation, albeit bond percolation. 

The return probability of a random walker can be calculated numerically with a simple com- 
puter program. One can let N w walkers perform at maximum N s steps, and count the number 
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of walkers that have re-visited the site they started from. Alternatively, one can let N w walkers 
perform steps until they are either back at the origin or further away from the origin than a 
certain predefined distance. I have used both methods to estimate P r with an estimated inaccu- 
racy of about 0.001. Note that there are two sources of inaccuracy. The first one is a statistical 
uncertainty, which scales as 1 / V N w . The second one is a systematic error, because each walker 
is stopped at a certain moment (after N s steps or at a given distance from the origin). For each 
of these walkers there is a finite probability that they would have reached the origin, when given 
enough time. As a result, the numerical estimates have a systematic error to the downside. This 
bias can be made smaller by using a large number of steps, or a large cut-off distance. In the 
calculation of the numbers quoted in Table [j], I used N s — 10 6 and higher, and N w — 10 5 . The 
results for the cubic lattices are consistent with Finch (1997) and Flajolet (1995). 

The numerical values for P r are close to the percolation thresholds, especially for d > 4. It looks 
as if the return probability of a random walker obeys the same scaling relation as the percolation 
threshold. In the remainder of this section I therefore calculate a crude approximation for the 
return probability using simple arguments. This approximation corroborates the scaling of the 
return probability as P r ~ l/d. 

Consider a random walker on the (i-dimensional Kagome lattice. Since all directions of the 
lattice are equivalent, it does not make a difference which step the walker makes first. Assume, 
without loss of generality, that the walker remains within the base polytope {i, 0}. Assume further, 
as a first approximation, that he stays within this polytope for a number of steps, and then returns 
to the origin. The walker can take an arbitrary number of steps within this base polytope, with 
a probability (d — l)/(2d) (the walker can choose from 2d directions, d of which are outside the 
base polytope, and one of the 2d is the origin). After a number of steps the walker should step 
back to the origin, which happens with a probability l/(2d). This approximation yields 

^ i (d-iy- 2 i , . 

^ = £^Ur) =TTd- (10) 

s—2 v / 

This approximation already reveals an interesting point: in high dimensions, the return probability 
scales as P ~ l/d. Since the probability 1/(1 + d) is the exact probability of return via a few 
selected paths, we also know that the exact probability on return via any path will be higher 
than 1/(1 + d). It is therefore impossible that it is as low as l/(2d — 1), as is the case for cubic 
lattices. 

One can improve on the above approximation by allowing the walker to step outside the base 
polytope occasionally. Choose for instance i sites x\ from which the walker steps outside. If the 
walker makes s steps within the base polytope before returning to the origin, there are ( s 7 M 
possible choices (because the last of the s steps is to the origin): 




The factor I /(id) emerges because the walker has a probability of 1/2 to step outside the base 
polytope, and a probability 1 / (2c?) to immediately step back to Xi. One can also allow the walker 
to make a number of additional steps, as long as he is only one step away from Xi and outside the 
base polytope. There are (d — 1) possibilities out of a total of 2d to make such a step. 

l l + 3/(2rf) 

d ' 1 + 3/ '(2d) + 3/ (2d 2 )' 

This type of reasoning can be taken a step further, yielding 

= 1 l + 2/rf + 9/(4d 2 

r ' 3 " d' l + 2/d+ll/(4d 2 ) + 9 
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All the paths that are included in this way are paths where the walker returns to the origin through 
the base polytope he started in. The first contributions from paths that return from the opposite 
side are paths of six steps, see Fig. ||. As there are 2d(d — 1) of such paths, the probability to 
return from the opposite side in six steps is (d — l)/(2<i) 5 , which is of the order 1/d 4 . Thus the 
conclusion that the return probability scales as 1/d remains unchanged. 

As a numerical check on this scaling behaviour, I calculated the return probability for the 8- 
dimensional Kagome lattice: P r = 0.126±0.001. This compares well with the approximation 0.124 
from Eq. Q. 



6 Discussion 



The Kagome lattice has rather different properties in d dimensions than the cubic lattice: it does 
not behave as a Bethe lattice in high dimensions, although one usually assumes that all lattices 
do. The series expansion of Section || provides some insight into why this is the case. One clear 
difference with the cubic lattice is that the leading power of q in the expression for the number 
of clusters of s sites is (s + l)d — (s — 1) for the Kagome lattice, but 2sd — 2(s — 1) for the cubic 
lattice. The coefficient of sd differs by a factor 2 here. The significance of this power of q is that it 
is the 'coordination number' of a cluster of s sites, just as q is the coordination number of a single 
site. The underlying reason for the different powers of q is that two adjacent sites on the Kagome 
lattice always have d— 1 common neighbours, whereas two adjacent sites on the cubic lattice have 
none. In high dimensions this difference does not disappear, but instead becomes more important. 

The same difference also explains why the return probability is different for these lattices. 
Since adjacent sites have d — 1 common neighbours, there is, in high dimensions, an increasing 
number of ways to return to the origin in three steps. On a Bethe lattice, such paths do not exist: 
here a walker can only return to the origin by retracing his steps. Therefore the Kagome lattice 
does not behave as a Bethe lattice in high dimensions. 

Going back to the 2-dimensional Kagome lattice, one could say that one of the reasons why 
it has a much higher site percolation threshold than the square lattice, is because of its 'common 
neighbours'. 

We can apply the same reasoning to bond percolation. The bond percolation threshold of the 
Kagome lattice is also higher than that of the square lattice, but the difference is smaller than for 
the site percolation threshold. Two adjacent bonds on the square lattice always have 2 common 
neighbours. On the Kagome lattice, two adjacent bonds have, depending on their relative position, 
2 or 3 common neighbours. The former happens in | of the situations, the latter in |. Therefore, 
for bond percolation the Kagome lattice has more resemblance to the square lattice, but still the 
two lattices are different. This is consistent with the difference in bond percolation thresholds 
being smaller, but not negligible. 

Finally, let us consider the relation to the diamond lattice in d dimensions. Since the site 
percolation threshold of the Kagome lattice scales as p c ~ 1/d = 2/q, this holds for the bond 
percolation threshold of the diamond lattice as well. But for the diamond lattice q = d + 1, so 
Pc,dia ~ — 1), which is the familiar Bethe result once more. This highlights a more general 
property. Each bond problem can be mapped onto a site problem (see e.g. Kesten 1982). When the 
bond problem is on a lattice with q = qj,, the lattice for the corresponding site problem will have 
q s = 2(q b — l). As a consequence, for each class of lattices for which the bond percolation thresholds 
scale as p c ~ /(<?), there is another class of lattices for which the site percolation thresholds scale 
as p c ~ /(|<7 + 1)- It would therefore be more consistent to describe bond percolation in terms of 
the number of neighbours a bond has, which is 2(q — 1). 
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a b 

Figure 1: (a) A description of the Kagome lattice as a lattice with a 3-point basis. The basis 
points form a (regular) triangle, (b) Three base triangles define a larger triangle, which truncates 
to a hexagon. 




Figure 2: The star-tetrahedron transformation. The dashed lines AX,. . .,DX, that come together 
in the point X, form the diamond lattice. 
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Figure 3: The scaling of the percolation threshold with network size N = (d + l)L d . 
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Figure 4: The percolation thresholds as a function of dimension d. Eq. (Q) is the result of Gaunt 
et al. (1976) for the thresholds of the cubic lattices. The thresholds for the Kagome lattices scale 
as 1/d. 
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Figure 5: The clusters of 3 sites. 
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Figure 6: The clusters of 4 sites. 




Figure 7: The clusters of 5 sites. 




neighbours on two sides of the cluster 



Figure 8: On the Kagome lattice the first instance of a cluster with fewer than usual neighbours 
occurs for clusters of five sites (left picture) . On a cubic lattice it occurs for clusters of three sites 
(right picture). 
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